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Abstract 

This paper presents the stabiUty analysis of a system sUding at low velocities (< 100/um.s~^) 
under a periodically modulated normal load, preserving interfacial contact. Experiments clearly 
evidence that normal vibrations generally stabilize the system against stick-slip oscillations, at least 
for a modulation frequency much larger than the stick-slip one. The mechanical model of Bureau 
et al. (2000), validated on the steady-state response of the system, is used to map its stability 
diagram. The model takes explicitly into account the finite shear stiffness of the load-bearing 
asperities, in addition to a classical state- and rate-dependent friction force. The numerical results 
are in excellent quantitative agreement with the experimental data obtained from a multicontact 
frictional system between glassy polymer materials. Simulations at larger amplitude of modulation 
(typically 20% of the mean normal load) suggest that the non-linear coupling between normal and 
sliding motion could have a destabilizing effect in restricted regions of the parameter space. 
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I. INTRODUCTION 



It is well known that non-linearities in the constitutive laws of dry friction may lead 
to the instability of steady frictional sliding against stick-slip oscillations, even for a single 
degree-of-freedom system driven at nominal constant velocity through a compliant stage. 
Sliding instability is an important issue in mechanical engineering since it is an ultimate 
limitation to the positioning accuracy for precision structures and machine tools. Moreover, 
stick-slip oscillations may be strongly non-linear and make servo-control extremely difficult. 
When designing a sliding mechanism, it is therefore of primary importance to characterize 
accurately the variations of the friction coefficient with, e.g., sliding velocity, keeping in 
mind that even slight variations may have a destabilizing effect. This requires to go beyond 
Amontons- Coulomb's law which assumes a constant friction coefficient. It might be legiti- 
mately feared that a more detailed constitutive law would have a restricted scope, e.g., in 
terms of materials and range of sliding velocities. It is thus remarkable that in the limit of 
low velocities (typically lower that lOO/xm.s"^), and light enough loads so that the inter- 
face is made of a sparse set of micro-contacts between load-bearing asperities, a relatively 
material-independent frictional behaviour is found which can be accounted for by a simple 
set of non-linear constitutive equations. Such studies have been initiated in the field of rock 
mechanics by Dieterich |I|] and Rice and Ruina [0], motivated by the need for low velocity 
friction models to investigate fault dynamics and earthquake nucleation. They have put on 
a firm phenomenological basis the idea, already suggested by the work of Rabinowicz [Q, 
that friction does not depend only on the instantaneous sliding velocity v but also on the 
whole sliding history. An experimental signature is the hysteretic frictional response of the 
interface when the slider is driven at a non-steady rate. Rice and Ruina proposed a family 
of dynamical equations coupling the sliding velocity to a set of state variables. Subsequent 
experimental investigations have shown that a single state variable is sufficient for most 
purposes. These experimental studies were performed on a wide range of materials, such 
as granite paper [Q, polymer glasses and elastomers 0. The friction force in this 
model is F = W fi{v, 0) with W the normal load and fi the friction coefficient. Moreover, 
it has been possible ^, to give a physical interpretation of as the average "age" of the 
microcontacts which grow while the material creeps under normal load, until sliding inter- 
rupts the process by renewing the load-bearing contact population. The dynamical model is 
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closed by specifying a differential equation coupling f to so as to account for the renewal 
of the microcontact population after a slip length Dq of micrometric order. This length is of 
order the mean radius of the microcontacts between surfaces of micrometric roughness 0. 
The resulting state- and rate-dependent friction laws will be hereafter referred to as SRF. 
Among several SRF expressions proposed originally, the one that we use in this paper are: 



lj.{v, 0) = /io + A In (^^) +B\n(^l + ^) (1) 
for the friction coefficient and 

for the evolution of the state variable, where /zq, A, B, Vsat and Vq are constants. 

This SFR model has been extensively validated by testing against numerous experimental 
situations involving transient dynamical responses of the system. The most stringent test 
relies upon the non-linear characteristics of the bifurcation from steady-sliding to stick- 



slip oscillations |Ty]. The model can be understood as resulting from two distinct physical 
mechanisms, the effect of which can be summarized in the following decomposition of the 



friction force, proposed by Bowden and Tabor |]TT|, in terms of the real area of contact J2r, 
and an interfacial shear strength as'- 



Fiv,^) = as{v)K{^) . (3) 

Here, the real area of contact depends on the interfacial age because it grows due to the 
creep of the load-bearing asperities . The velocity dependent interfacial shear strength has 
been ascribed to the adhesive, nanometer thick junctions between microasperities. A simple 
microscopic model has been proposed for the elasto-visco-plastic rheology of the junctions, 
compatible with the existence of a finite friction threshold [§. 

Recently, attention has been paid to the effect of a time-dependent normal load on the 
response of a single degree-of-freedom sliding system. This situation is of practical inter- 
est when the mechanical design allows cross-talking between the normal (loading) and the 
tangential (driving) forces [12], or when external vibrations contribute to the loading of the 



interface, as it may be the case for seismic faults |13, 14, 15 1. The response of the system 



is not intuitive. First, since the friction force is directly proportional to the normal load, 
the sliding velocity is dynamically coupled to the normal load modulation, hence feeds back 
the friction force. Moreover, it has been shown that more subtle interplays must be taken 
into account. Linker and Dieterich []T3[ have interpreted the transient response to a step 
in normal load by coupling directly the time-variations of and W, thus adding a term 
— const (f)d\n{W)/dt in Eq. (^. The physical motivation for this extension of the SRF 
ageing equation is the fact that, according to P], a change in normal force creates fresh 
load-bearing contact area. This certainly influences directly the age 0, though probably 
in a weaker measure than proposed by Linker and Dieterich as briefly discussed in 
T6| . More recently. Bureau et al. [|16]| have studied the response of a sliding system to a 



periodically modulated normal load W(t) = Wq[1 + ecos(ti;t)] with e < 1. They found that 
the friction force, averaged on a modulation period, is significantly lowered with respect to 
the situation under constant load Wq. The oscillating part of the force, primarily reduced 
to its spectral component at u in the limit of vanishing e, becomes quickly anharmonic as 
e is increased while still remaining much smaller than 1. They have shown that the SRF 
equations can fit accurately all their results provided that the model is modified to ac- 
count for the finite interfacial shear stiffness k resulting from the elastic deformation of the 
load-bearing asperities. This means that the sliding velocity differs from the velocity of the 
center of mass of the slider, a statement which is clearly illustrated in the static state, i.e., 
for tangential forces well below the static threshold, where the interface responds elastically 



without sliding [1^, Under constant normal load and constant driving velocity V, this 
"hidden" interfacial degree of freedom manifests itself only for non-steady motion, and plays 
no significant role at the circular frequency Qc ~ V/ Dq of the oscillations at the onset of the 
stick-slip instability. 

However, under a modulated load at u; 3> one must take the finite interfacial com- 
pliance K into account, all the more so since the latter is known to be itself proportional to 
W [|l^. This results in a non-trivial and efficient coupling between the normal load and the 
sliding velocity. 

Of particular interest is the effect of load modulation on the sliding stability of the system. 
Dupont and Bapna [|12] have computed the critical stiffness of the drive below which stick- 
slip occurs for a slider-spring loaded at a constant angle with respect to the sliding plane. 
This configuration would provide a direct test for the coupling between (p and W proposed 
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by Linker and Dieterich |jT3[, but the experimental study has not been performed so far. 

The present paper addresses the problem of the stability of a slider-spring system under 
an externally and harmonically modulated normal load. The experimental arrangement is 
described in section 2 and it is shown that for a circular frequency uj ^ the modulation 
generally stabilizes the system against stick-slip. This spectacular effect is accounted for by 
the SRF model with modulated interfacial stiffness, as shown by the numerical study of the 
bifurcation which is detailed in section 3. 



II. EXPERIMENTS 



A. Apparatus 



The apparatus (Fig. |l|) consists of a slider of mass M driven along a track through a 
loading spring of stiffness K = 0.21 N./im^^. The loading end of the spring is moved at 
constant velocity V in the range 1-100 /im.s~^ by means of a translation stage driven by a 
stepping motor. The spring elongation is measured by an inductive probe (Electro, sensor 
4937, module PBA200), with a 0.1 fim resolution over the 10 kHz bandwidth. The average 
normal load Wq can be set in the range 3-23 N with the help of a vertical spring attached 
to a remote point itself translated horizontally at the pulling velocity V through a second 
translation stage, in order to prevent any tangential coupling. The normal load modulation 
is achieved by means of a vibration exciter (LDS, model VlOO) rigidly attached to the slider: 
a harmonic voltage input of given amplitude and frequency / results in a harmonic vertical 
motion of the moving element of the exciter on which an accelerometer (Briiel & Kj^r, 
type 4371 V) is fixed. An acceleration of amplitude 7 of this moving element of mass m 
induces a normal load modulation on the slider of amplitude rwy at frequency /. We thus 
obtain a normal load W{t) = Wo[l + ecos{ujt)] with uj = 27r/ and e = mj/Wo in the 
range 0.01-0.5. A fixed frequency / = 120 Hz has been used for the whole study. Two 
poly(methylmethacrylate) (PMMA) samples are glued, respectively, on the slider and the 
track. They have nominally flat surfaces which have been lapped together with 400-grit 
SiC powder and water to obtain a rms roughness Rg = 1.3 fim [|l^. The interface between 
the two blocks is made of a sparse set of load-bearing microcontacts ll8|. An air layer of 



micrometric thickness is therefore trapped between the surfaces and acts as a viscoelastic 
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element, in parallel with the microcontacts, which partially bears the normal load. This 
effect has been studied in details in Bureau et al. |T6[ who concluded that the remaining 
effect of the load modulation on the asperities can be described by an effective amplitude 
^eflf = P^i with p a constant close to 0.5, taken in the following as p = 0.43, a value which 
will be justified in section [111 C] . 

B. Localization of the stick-slip bifurcation, effect of the modulation 

The bifurcation between steady sliding and stick-slip oscillations under constant load 
(e = 0) has been extensively described (see e.g. [H). When K and V are kept constant, 
steady sliding occurs for values of the remaining control parameter Wq < Wq{V) where 
the i^'-dependency has been omitted here since the value of K is fixed in this study. The 
bifurcation is of the direct Hopf kind, which means that the amplitude of oscillation of the 
slider velocity is vanishing ^ when approaching Wq from below {i.e. increasing Wo up to 
Wq), while the circular frequency tends to a finite value Qc- In addition, the characteristic 
time of the oscillating transients diverges when approaching Wq from above {i.e. decreasing 
Wq towards Wq). A practical consequence is that as the bifurcation is tracked down, it 
becomes increasingly difficult to distinguish between steady stick-slip oscillations and tran- 
sient relaxation towards steady sliding, resulting from the perturbating effect of friction force 
fluctuations along the track. 

For e = 0, the ratio K/Wq is the relevant control parameter, at least in the low velocity 
region where inertia of the slider oscillating at the circular frequency Qc can be neglected . 
Henceforth, although the external stiffness K is kept constant for the whole set of experi- 
ments reported in this paper, we will keep on representing the stability domain of the system 
in the parameter plane {K/Wq, V) where it is bounded by the experimental bifurcation curve 
(Fig. 1^). The experimental uncertainty on the critical value K/Wq, determined from the 
standard deviation over at least 10 measurements at a given velocity V, is indicated by the 
error bars in Fig. 0. It is typically ±3%, except at the larger velocity, since the results are 
more sensitive to long wavelength irregularities along the track for large sliding distances. 

When a harmonic modulation at u is superimposed to a value of Wq corresponding to 

^ Note that the term "shck-shp" is therefore a misnomer since the shding velocity does not reach zero, i.e., 
the shder does not "stick" during an oscillation period. 
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steady sliding at V, the velocity of the slider oscillates about V, possibly in a anharmonic 
way, with a fundamental component at uj. The motion of the slider, when averaged over a 
period 2tt/uj is therefore steady. For given V and K, one has now to consider two control 
parameters, namely Wq and eefr- 

When the d.c. load Wq is increased while keeping constant e^s, the average slider motion 
is found to become of the stick-slip kind above a critical value, VFJ (V^, Ceff) > W^o (^) (Fig- 0)- 
A normal load modulation of even very small effective amplitude may therefore stabilize 
the system against stick-slip as illustrated directly in Fig. || where a modulation with eefr = 
4.5 X 10~^ is enough to suppress well developed, strongly anharmonic, large amplitude and 
low frequency stick-slip oscillations (the force signal then only shows the the remaining small 
amplitude modulation at the forcing higher frequency). 

The effect of ecs on the critical value of K/Wq{V) is characterized in Fig. ^ The higher 
the velocity, the stronger the stabilizing effect of the normal load modulation. The effect is 
spectacular when described in terms of the velocity domain corresponding to steady sliding 
at constant K/Wq. For instance, the critical velocity at K/Wq = 0.026 /xm~^ is decreased 
by more than a factor of ten by applying a modulation with tcs = 0.09. 

The empirical study indicates that, as a rule of thumb, steady-sliding is promoted by 
high velocity, high amplitude of normal load modulation, low average normal load and large 
stiffness. This is tested in the following against a numerical study of the SRF model including 
normal load modulation. 



III. NUMERICAL STUDY 

A. The state- and rate- dependent fiction (SRF) model equations 

The SRF laws (Eqs. [1|, ^ are incorporated into the equation of motion of the slider, 
according to the simple model sketched in Fig. ^. The proportionality constant between the 
normal load W and the interfacial stiffness k is a length A of micrometric order: k, = W/ A 
Tsj . The equation of motion of the slider thus reads: 



Mx" = K{Vt -x)- k{x - Xpi) (4) 



7 



where here and henceforth the prime denotes time derivative. Taking a massless interfacial 
zone (a reasonable assumption, see appendix), we also have, according to Eq. (|I[): 



W 
T 



(x — Xpi) = W 



+ 5 In 1 



sat ' 



(5) 



where v = x'^i is the relative sliding velocity at the interface and follows the evolution law 
(§) rewritten here for the sake of clarity ^ : 



Do 



(6) 



For numerical purposes, we wish to recast those equations in the form of a system of 
first order ordinary differential equations (ODEs). Noting z = x',u = Vt — x, further 
differentiating x — Xpi with respect to time in (H) using the exphcit expression for W(t) = 
Wo[l + ecscos{ujt)] and solving for v', we get the following ODE system, which we will use 
for the numerical bifurcation analysis: 



u 



V-z 
K 



M 



:U 



Wo{l + ecscos{ujt)) 



V 

1 - 



M 

z — V — 

V (j) 



/Wo + ^ In I — 



AEKat(l-^0/I^o) 
/^O(l + Kat0//^O) . 



5 In 1 + 



sat ' 







(7) 

(8) 

(9) 
(10) 



B. Determination of the SRF parameters 

In order to analyse the data within the SRF framework, we need to determine a set of 
values for the relevant parameters of the model. This is performed under constant normal 
load, according to a well established procedure. The values, which will be used in the 
numerical analysis, some of them as trial ones, are gathered in Table |I[ The useful formulae 
are established in the appendix. 

^ This equation, including the actual sliding velocity v in place of x' is consistent with our physical under- 
standing of the age 4>. It has been checked that mistaking x' for v, as in jT6|, has no significant effect on 
their results, at least for the modulation frequencies much larger than V/Dq used in their study. 
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(i) First, the steady sliding friction coefficient fidiV) is measured to be velocity-weakening 
with an almost constant logarithmic slope over the 1 - 100 /im.s"^ range. This indicates 
that Vsa.t, above which ^diV) increases with increasing V according to Eqs. (|1], 0), is 
certainly larger than 100 /im.s^^ and allows to extract a value ioi B — A and fiQ at 
Vq = 1 /im.s^^, i.e. far below the saturation of the ageing term. 

(ii) Next, the critical value of K/Wq, here taken at midrange (10 /im.s^^) where inertial 
terms can be neglected in Eq. (pTD, yields a determination of the memory length Dq. 



(iii) From the value of the critical stick-shp period at 10 /xm.s ^ (Eq. (0)), we obtain a 
value for A. 

(iv) A determination of Kat is finally obtained by a best fit of the whole bifurcation curve 
in the plane [K/Wq, V) for V in the 1-50 fim.s~^ range, treating the inertial term 
in Eq. (^) as a perturbative one. Since the value of V"sat is out of the experimental 
velocity window, this determination is not very accurate. Treating several data set 
corresponding to different runs yields an uncertainty as large as ±25% on the value of 

(iv) The value of the length A, defined by the ratio of the load Wq and the interfacial shear 
stiffness k, has been obtained in [TBI from a best fit of the a.c. response of the slider 



position to the normal load modulation. 

It is clear that this procedure, though systematic, generates cumulative errors which are 
difficult to evaluate (the uncertainties on A, B and Dq given in Table (|) are conservative 
values). In view of the high sensitivity of the bifurcation to small variations of the parame- 
ters, we have chosen in the following numerical analysis to use the set determined above as 
a trial one. Namely, the parameters which are left free are A, i?, Dq, A, Vsati and the ratio 
p = eeflf/e. 



C. Bifurcation analysis 

Technically, the transition from steady sliding to stick-slip, both states being modulated 
by the forcing (when e 7^ 0), is a Neimark-Saker bifurcation (also called secondary Hopf), 
which corresponds to two complex conjugate values of the fundamental matrix of the ODE 
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system crossing the unit circle. The fundameltal matrix H is defined as dH/dt = J{x{t))H 
with H{x, 0) = I, with J the Jacobian matrix of the ODE system and I the identity. The 
numerical software CANDYS/QA |]19[ has been used to track this bifurcation. For a given 
parameter set a bifurcation curve like in Fig. ^ can be obtained as follows: for a given 
driving velocity, one starts from a low enough normal load Wq in order to be in the steady 
sliding regime. Once such a "first point" is indeed found by the software, one varies Wq only 
(1-parameter continuation) until a bifurcation is detected {Wo = Wq); one then follows this 
bifurcation curve by further varying the driving velocity too (2-parameter continuation). 
The procedure to detect the bifurcation has been automated. Starting from the param- 



eter values determined in |III B| , the critical values Wq are determined and compared to the 
experimental ones Wq^^^. A systematic procedure (Powell's method, as described in p^) is 
then used that attemps to minimize '^vel^oiP^) ~ ^Oexp]^' with p = eefr/e, for a represen- 
tative set of experimental data. The set of parameters hence determined, which corresponds 
to a local minimum of the cost function, reducing it by a factor of fifteen ^, is given in 
Table |. 

The full bifurcation curves are then determined as described above, by a 2-parameter 
continuation. The results are shown in figures ^ and ^ together with experimental data 
covering a range of ees values wider than the one involved in the adjustment procedure. 



IV. DISCUSSION 

The curves displayed in Figs. ^ and ^ have been calculated with the optimized set of 
parameters. However, it is worth noting that the trial set yields numerical results in good 
qualitative agreement with the experimental data as well. Namely, the main effect of normal 
load modulation, at least for moderate values of eefr, which is to stabilize sliding against stick- 
slip oscillations, is well reproduced by the SRF model. Moreover, the enhanced efficiency of 
the modulation on increasing the sliding velocity is correctly accounted for. 

The set of optimized data differs from the trial set essentially for three parameters, namely 

^ When several critical values for Wq can be detected, the one retained in the evaluation of the cost function 
is the closest chosen among the odd ones (first, third, etc.), corresponding to a transition from steady 
sliding to stick-slip when increasing Wq, while even ones correspond to a transition from steading sliding 
to stick-slip. 
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the elastic length A, the saturation velocity V^at, and the ratio p = ecg/e accounting for the 
air-cushion effect. 



The final value for A lies within the error bars estimated in HTB 



As already mentioned in Sec. [Ill lj| , the large variation of Kjat during the optimization 
procedure is attributable to the fact that the crossover from a velocity- weakening regime to a 
strengthening one for steady sliding friction lies well above the upper experimental velocity, 
whence the goodness of the fit is only weakly sensitive to V^at- 



The ratio p = eefr/e was determined in |16| by comparison of the experimental shift of the 
steady friction level Apo at 120 Hz and the value predicted by the SRF model. The value 
taken in this reference was p = 0.48. Taking into account the error bars on A/iq one finds 
that the relative uncertainty on p is about ±20%. The optimized value for this parameter 
lies therefore within this range. 

Thus, the SRF model with its set of parameters as determined from the dynamical study 
of the system under constant normal load is fully predictive as regards the sliding stability 
of the system under modulated load, at least for the values of ees probed by the data of 
Fig. 0. In turn, the sensitivity of this experiment enabled us to refine the determination of 
the parameters. 

The quantitative overall agreement between the experimental data and the numerical 
curves in Fig. ^ is excellent for, say, e^s < 0.15. Above this value the calculated curves 
tend to fold and correspond to a re-entrent stability diagram; namely, for given ees and V, 
increasing K/Wo yields successive bifurcations from stick-slip to steady sliding then back 
to steady-sliding, etc. No experimental evidence of such an unexpected behaviour has been 
encountered so far. It is clearly the result of the non-linear coupling between the normal 
load modulation and the stick-slip oscillations. As such, it is expected to depend drastically 
on the details of the SRF laws. The importance of the terms which ultimately cut-off the 
logarithmic variations in the SRF laws has been stressed in several studies [|10], The 
existence of V^at, which accounts for a short time cut-off in the creep deformation of the 
load-bearing asperities 0, yields one of these terms. It should be kept in mind that the 
SRF constitutive law ([1|) retains only the leading terms in the expansion of the friction 
force in powers of ln{v). For instance, Eq. (|^) with physically sounded expressions for 
(Ts{v) and would lead to terms of order AB ln{v) ln{(j)) which, though negligible for 

most purposes, would probably affect the critical behaviour of the system under a strongly 
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modulated load. For these reasons, we think that a full quantitative agreement between the 
experimental bifurcation at large ^ ees and the SRF model predictions would be illusive. 

In addition, we have investigated numerically the effect of the extra term in Eq. 
proposed by Linker and Dieterich []T3[ and already discussed in [|T^. No significant effect 
has been found at 120 Hz, where we conclude that this term, if it exists as such, is not 
relevant to the present relatively high frequency study. 



V. CONCLUSION 



The stability of a sliding system with a few degrees of freedom, submitted to a periodi- 
cally modulated normal load, has been studied experimentally. The study clearly evidences 
the role of load modulation, even at moderate amplitude, as a stabilizer against stick-slip 
oscillations. Though this effect of vibrations is seemingly part of the empirical culture in 
mechanical engineering, it is the first time, as far as we know, that it is investigated experi- 
mentally. This effect should be of great interest in fault mechanics as well as in the control 
of precision structures. The results have been compared to the numerical predictions of a 
model of the SRF type, relevant to multicontact friction at low velocities and low loads, 
including finite interfacial shear stiffness as a key parameter. Excellent quantitative agree- 
ment has been found as long as the amplitude of load modulation is restricted to about 10% 
of the dead load. 

Although, as discussed above, the main effect of the normal load modulation is stabiliza- 
tion, the numerical study strongly suggests that destabilization may also occur, due to the 
highly non-linear features of the model which also gives rise to rentrent stability diagram 
in Fig. 0. More precisely, it can be seen in this figure that the = 0.18 curve crosses the 
eeff = 0.13 one around V = 7/im.s~^. For a {V,K/Wo) point slightly on the right of this 
crossing, in between the two curves, increasing tcs would result in a bifurcation from stable 
sliding to stick-slip. This effect has not been observed directly so far, probably because it 
corresponds to small regions of the parameter space, strongly dependent on the value of the 
parameters. Clearly, this point would deserve further experimental study. 

^ It has been shown in [ p6[ that the relevant perturbation parameter is actually fj-o^osf-A ^ Eoff which is 
already larger than 1 for Ccff = 0.1. 
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LINEAR STABILITY ANALYSIS FOR e = 



The linear stability analysis of the SRF equations has been performed previously p], 
However, difference between sliding velocity and the velocity of the center of mass of 
the slider was disregarded in these works. Since the interfacial stiffness is of paramount 
importance when the normal load is modulated at relatively high frequency, it is necessary 
to evaluate its role on the location of the bifurcation under constant load. Moreover, we will 
derive in this appendix expressions for the critical stiffness and the critical pulsation that 
hold for any state- and rate-dependent friction force. 

Let us consider a general expression: 

F = Wo^Ji{v,ct>) . (11) 

The time-evolution of the age variable is ruled by: 

<^' = l--f. (12) 

When the slider is driven at constant velocity V , the steady sliding values of the dynamical 
variables are v = V and (p = Dq/V. We define: 

The position of the center of mass of the slider is x{t) so that the elongation of the loading 
spring is a; — Vt. At frequencies of interest the interfacial zone can be assumed massless and 



essentially elastic with a frequency independent, real stiffness k [|T^. The following relation 
thus holds: 

, d f Mx" - K{yt-x)\ , 

^ = ^+11 -\ ■ (14) 



dt\ K 

We will make use in the following of the ratio r] of the loading spring stiffness K to the 
equivalent stiffness of the loading spring in parallel with the interface || k: 

K 

Finally, the dynamical equation for the motion of the slider reads: 

Mx" = K{Vt -x)- Wofi{v, (j)) . (16) 
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The set of dynamical equations ([T^JT^JT6| ) is closed and can be linearized about the steady 
sliding state, setting: 



x = Vt- F{V, Do/V)/K + Sx, \5x\ < F{V, Do/V)/K 
<P = Do/V + 6<P, M-^Do/V . 

The linearized system becomes: 



(17) 



MSx" = -K5x - Wo [{lijV) {7]6x' + fSx'") + {fi^^V/Do 



-r]6x'/V - §6x"'/V - 5(j)V/Do 



(18) 



The solutions are the real parts of the complex 6x = 6xoexp{iQt) and S(j) = S(f)oexp{iQt) 
with Q a complex number. Replacing into Eq. ([18D and writing the condition for non-trivial 
amplitudes 6x0 and 500, one finds the dispersion relationship 



C4fi^ + Csn^ + Cafi^ + Cifi + Co 







(19) 



with: 



Co 



KV 
Do 



KDq 
Wo 



c. 



2 V ^0 



(20) 



kDo 



The critical value of the control parameters and the critical pulsation Vic are obtained by 
expressing that at the Hopf bifurcation (at least) one root of Eq. (^) crosses the imaginary 
axis. Setting that Vt is purely imaginary and extracting the real and imaginary components 



from Eq. (|T9D yield the requested values. 

Let us first solve for an infinitely stiff interface, i.e. for r] = 1 and k, 
straightforward to find: 



00. It is then 



(21) 



and 



Dq^c 
V 



(22) 
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These relations make sense only for fifp — fiv > 0, i.e. when the steady sliding friction 
coefficient yU***(y) is velocity weakening: dfi'^^ /dlnV < 0. For the particular expression of 

fj,{v, 0) used in the numerical analysis, this reads: 

du^^ B 

IJ^-IJ,= ^ , ^^.^^ - A > . (23) 



d\nV 1 + 

Now, let us evaluate the contribution of the finite interfacial stiffness k to Eq. by 

estimating the order of magnitude of c = |C4r2^/Co| with VLc given by Eq. (^), i.e. VLc — 

V/Dq. This reads: 



where we have expressed that WQ/{nDQ) = X/Dq ~ 1. One can estimate c < 10^'^ within the 
experimental velocity range, hence the fourth order term in Eq. (p!9D can be safely discarded. 
Next, a finite n introduces perturbative terms in C3 which are of order fiy, /x^ ~ 10~^, 
still well below the relative uncertainty on the experimental determination of the critical 
parameters. Since the correction to the drive stiffness due to the interfacial elastic element 
is of order 77 — 1 = K/k, ~ 10^^, it can be concluded that for the purpose of calculating the 
values of the critical parameters, the finite interfacial stiffness has no practical effect, and 
one can make use of Eqs. ( ^l]) and (|22|) . 
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FIGURES 




FIG. 1: Main elements of the experimental setup: Translation stage (Drv); Loading leaf spring 
(Lsp); Displacement gauge (Gg); Vibration exciter (Vb); Weighting spring (Spr); Accelerometer 
(Acc). 
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FIG. 2: Stability diagram for different values of the modulation amplitude. For given V and ecff, 
bifurcation from stick-slip to stable sliding occurs when the control parameter K/Wq overcomes 
the plotted critical value: ees = 0(«); 0.045(A); 0.09(H); 0.13(0); 0.18(a). For the sake of 
clarity, typical standard deviations are plotted as error bars only for eeff = 0. The solid line curves 
are the output of the numerical study (see section IIIC| ). The larger eeff the lower the curve at 
V = l//m.s~^. 
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FIG. 3: Time evolution of the loading spring elongation for V = 8 nm.s~ and different modulation 
amplitudes eefr indicated at the right end of each trace. A vertical offset has been added to each 
trace in order to display clearly the bifurcation sequence from stick-slip to stable sliding. The inset 
is a blow-up of the stable sliding trace showing the remaining oscillating response at the frequency 
of the load modulation (/ = 120 Hz, much higher than the stick-slip frequency). 




FIG. 4: Reduced critical load vs. Ecs for different driving velocities: V{fim.s^^) = 1(a); 5(«); 
10 (□); 30 (■); 50 (o). The curves are the output of the numerical study (see section III Q ) , labeled 
with the corresponding velocities in fim.s"^. 
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FIG. 5: Equivalent mechanical circuit of the slider/track system. K is the stiffness of the loading 
spring, K, is the one of the interface. 



TABLES 



TABLE I: Values of the SRF parameters. 



parameter 


trial value 


optimized value 


K (N./xm-i) 


0.21 




M (kg) 


2.37 




/xo(Vb = 1 Atm.s"^) 


0.33 




A 


0.012 ± 0.002 


0.0126 


B 


0.023 ± 0.002 


0.0241 


Do (/^ni) 


0.40 ± 0.04 


0.402 


V;a,. (/au.s-1) 


280 ± 70 


25G 


A (/xm) 


0.62 ±0.15 


0.56 


p = eeff/e 


0.48 ±0.10 


0.43 
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